Nucleation model for the description of glass formation 
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We present in this letter a model of glass formation using energy barriers and a crystal nucleation 
process. We then analyze the corresponding dynamical equation in the vicinity of the stationary 
solutions. The occurence of a pure amorphous solution (i.e. glass) is due to the contribution of a 
ratio A depending on the cooling rate and the crystal nucleation frequency. We also construct time- 
temperature transformation (T-T-T) curves in order to relate the model with the kinetic treatment 
of glass formation. 
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The formation of glasses requires cooling to a suffi- 
ciently low temperature, below the glass transition tem- 
perature, without the occurence of detectable crystalliza- 
tion. 

One of the most comprehensive theory of dynamical 
processes in glasses is the so-called mode-coupling the- 
ory (MCT) , developed by Goetze and others 0. The 
MCT is based on a set of non-linear damped second- 
order differential equations, which couple the density cor- 
relation functions of a supercooled liquid in a retarded 
way. These equations have several interesting properties, 
such as the possibility of a stretched-exponential time 
decay for the correlation function for particular values of 
the strength of the non-linear coupling. They exhibiti, 
also, a kind of " ideal glass" transition, at an ideal glass 
temperature Tc, corresponding to a singularity of the 
equations. Nevertheless, the major limitation is that, so 
far, only spherically symmetric interatomic interactions 
have been considered, as for instance the computer cal- 
culation by Bengtzelius on a Lennard- Jones system [Q. 
Thus, the complexities in real glass-forming systems aris- 
ing from covalent interactions (e.g. silica) are neglected, 
so that the theory seems only well-adapted for fragile 
glass- forming liquids (e.g. organic or molecular glasses). 
There have been recently also theoretical developments 
of mean-field models in order to understand the role of 
different factors in the glass transition, such as config- 
urational entropy [gj, energy barriers |I| or free energy 
||, some of them being inspired from spin glass models. 
However, the critical question in discussing glass forma- 
tion is not whether the solid will be obtained from a melt 
quench, but how fast a given liquid must be cooled in or- 
der to prevent the kinetic processes involved in crystal- 
lization. Thus, the cooling rate and the crystal nucleation 
rate should play the major role. 

The work presented here consists in the construction of 
an exactly solvable mean-field model of crystalline nucle- 
ation, based on a probabilistic description, which can be 
described by a dynamical equation. Although this model 
is rather abstract and has no thermodynamic transition, 



we believe that it captures at least parts of the physics 
involved in the glass transition and explains it in a simple 
fashion. The model attempts also to relate factors that 
are widely viewed as decisive in the formation and prepa- 
ration of glasses, namely structural, kinetic and thermo- 
dynamic factors ||. We shall see that the probability of 
occurence of an amorphous phase (identified with a su- 
percooled liquid below the melting temperature and then 
with a glass when temperature has decreased enough) 
depends strongly on the cooling rate. Let us try to de- 
scribe what happens on a microscopic scale when a liquid 
is cooled and translate this in terms of thermodynamics. 

A liquid at high temperature (with T higher than the 
melting temperature) should contain a lot of aggregates 
of different sizes and shapes. The atoms which com- 
pose some of the aggregates may be already organized 
with periodic character (identified with a crystalline nu- 
cleus), whereas others do not exhibit this feature and may 
display a disordered structural organization (amorphous 
structure). 

Construction of the model. Let p c = p(Cq) be the prob- 
ability of finding an atom, in a crystalline aggregate, be- 
ing in the initial liquid configuration Co- The structure 
associated with this aggregate is size-independent and 
there is no reference to a spatial structure, therefore the 
model is of mean-field type. In this description, the prob- 
ability p c corresponds to the number of atoms which are 
already trapped inside a crystalline unit cell or a set of 
such cells, over the total number of atoms. Similarly, 
we define the probability of finding an atom inside an 
amorphous aggregate by 1 — p c . The glassy state will 
thus correspond to p c — at low temperature and the 
crystal to p c = 1. The probability p c depends at least 
on time, temperature and the cooling rate. When the 
temperature is decreasing, the aggregates start to grow 
by nucleation. The energy barriers related to such ag- 
glomeration processes depend, of course, on the nature 
of the clusters which are involved. Therefore, we define 
for a crystal-crystal nucleation an energy barrier of 2E C , 
for a crystal-amorphous structure agglomeration an en- 
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ergy barrier of E c + E a and for an amorphous-amorphous 
structure agglomeration 2E a . E c is the mean energy 
stored in a crystalline cluster and E a the one stored in an 
amorphous structure. This means that the initial config- 
uration Cq can be considered as a system with two energy 
states E c and E a and that the atoms can be treated as 
particles able to occupy one of these states. After ag- 
glomeration or nucleation, the system is in a final con- 
figuration C. It is composed of three energy states (2E C , 
E a + E c and 2E a ), and the repartition of the particles 
can be computed. The probability of finding an atom 
trapped in such a new state may be proportional to the 
product of the initial probabilities (i.e. to the number of 
related particles) and a Boltzmann weight, inspired from 
the Thouless- Anderson-Palmer (TAP) description |t]] of 
spin glass theory (e.g. the pure crystal state with prob- 
ability p cc is proportional to p 2 c ). Such a description has 
been also used in a growth model of polygons § and in 
structural glass models M. 
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where a = [3(E C — E a ) with j3 the reverse temperature. 
Since the probabilities are all normalized by the factor Z, 
only the difference a is relevant. The normalizing factor 
Z, which can be regarded as a TAP-like partition sum 
ensures that J2i j Pij = 1j as ^ should be. During nu- 
cleation, an atom trapped inside a crystalline state, can 
be found in the nucleated configuration C (energy states 
2E C or E c + E a ) or still in the initial configuration Cq 
(energy state E c ). Its probability is expressed by: 



Pc(t) = 



1 - s(t) 



p(C ) + s(t) p(C) 



where: 



P(C) 



2p c 



The function s(i) with < s(t) < 1, introduced above de- 
pends on time and represents the growth rate (i.e. a mea- 
sure of the number of nucleated crystalline aggregates 
which are already formed). Its derivative s(t) = ds/dt 
is then the nucleation frequency or growth speed. When 
s(t) = 1, the nucleation process is finished. 

The derivative of the probability p c {t) obeys the fol- 
lowing dynamical equation, which is constructed from 
the equations above and which takes into account the 
effect of the cooling rate Q = —dT/dt. 



dp c (t) 
dt 



P(C) - P(C ) 



Q 



d P (C) 
dT 



At finite temperature and for long enough times, the 
model reaches a stationary state, corresponding to ther- 
mal equilibrium, thus we shall be interested in the sta- 
tionary solutions of the dynamical equation. We can 



solve this equation in terms of the probability p(Cq), since 
the probability p(C) of the final state is also constructed 
with them. This should happen in the liquid and super- 
cooled state for very slow cooling rates and in the final 
solid for T ~ 0. Obvious solutions of the equation are 
p{Co) = (pure amorphous state) and p(Cq) = 1 (pure 
crystal state). Besides solutions without physical mean- 
ing (p(Co) > l)j an intermediate solution exists which 
shifts with the cooling rate Q, the crystal nucleation fre- 
quency s(t) and the temperature: 
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where A = Q/s(t) and u — a/2. The existence of this in- 
termediate solution depends on A, which must be inside 
region II of fig. 1. 

Dynamics of crystal and glass formation. In order 
to see in which direction and under what condition the 
system can evolve, we have performed the linearization 
of the dynamical equation in the vicinity of the station- 
ary solutions and computed the corresponding relaxation 
times tq(T) and Ti(T), which characterize the conver- 
gence of the nucleation process towards equilibrium. The 
behavior of the liquid is mainly driven by the factor A 
and typical situations occur, according the sign of a. 




FIG. 1. A plot of the critical factors |Ai| (solid line), which 
yields pl nt = 1, and |A 2 | (dashed line), which yields p™* = 0, 
as a function of temperature T with E c — E a = —1. The 
intermediate solution p 1 " 1 exists only if |Ai| < |A| < |A2|. 
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When a < (the crystalline energy barrier being 
lower than the amorphous one), there are only two so- 
lutions and the p(Co) = 1 one is an attractor. This 
means that every fluctuation £ in the probability which 
could occur during the nuclcation process, tends to van- 
ish with increasing time. The dynamical equation be- 
comes: Ti(T)d£ = — £dt with relaxation time: 



a process at thermal equilibrium (i.e. when A is negligi- 



n(T) 



tanh - 
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aA 
T 



+ 4) tanh § + 



aA 
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Growth is then made only out of atoms being part of 
crystalline aggregates. On the contrary, it is easy to ver- 
ify that a fluctuation 77 in the vicinity of the amorphous 
solution p(Cq) = will grow exponentially with increas- 
ing time and the relaxation time tq(T) [To{T)drj = rjdt]: 



ro(T) 



1 + tanh : 



(«A _ 4) tanh a _ aA 



So, there is no dissipation of an amorphous fluctuation 
77. We believe that this situation is typical of systems for 
which there is no possibility in glass formation (e.g. all 
the elements except selenium, sulfur and tellurium). 




FIG. 2. The relaxation time to(T) as a function of temper- 
ature for two different processes (cooling [circles] and heating 
[circles and solid line]) and two different values of |A| (5, 10). 
The solid line corresponds to T equ . 



An interesting feature of this dynamical analysis is 
given by the relaxation times t (T) and T\{T) which ex- 
hibit a very different behavior following the process (cool- 
ing or heating). We define by T equ the relaxation time of 
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1 — e a . For low temperatures (T — > 0), 



the dynamics of the process is close to equilibrium since 
Ti(T) ~ r equ . Starting from a high T liquid, one observes 
that Ti decreases when temperature goes down. This 
means that it becomes more and more easier for the sys- 
tem to fall on the crystalline attractor. t\ exhibits a min- 
imum at a certain value T* = -A + yj A 2 - A(E C — E a ), 
which we identify with the melting temperature. It is the 
temperature at which the driving force towards the crys- 
talline aggregation (p(Cq) = 1) is at its maximum. In the 
heating process (A < 0), the p(Co) = and p(Co) = 1 
solutions are both attractive, the latter with a relax- 
ation time T\ ~ T equ for low T, but this solution diverges 
rapidly in the vicinity of a temperature Ti close to T*. 
For T > T\, the crystal solution is unstable and the amor- 
phous (liquid) solution only is attrative. In the range 
< T < Ti , the intermediate solution exists but it is un- 
stable. The described dynamical behavior shows that a 
crystalline solid which is heated from a low temperature, 
remains the preferential structure up to a temperature 
Ti close to the temperature T*. There, the crystal relax- 
ation time Ti diverges, which allows the possibility for the 
system to investigate other possible structural pathways, 
such as the amorphous one. 

When a > 0, there is a possibility of formation of an 
amorphous solid, because thep(Co) = solution is attrac- 
tive, whereas, the crystalline state is repulsive. In this 
situation, tq(T) displays the following behavior: starting 
from a high T liquid phase, one can easily check that 
the fluctuations in the vicinity of the crystalline solution 
grow exponentially with a relaxation time T\(T). The 
To(T) relaxation time is negative and close to 0. It starts 
to decrease when the temperature decreases (fig. 2, cir- 
cles), i.e. the dynamics of the attractor becomes slower, 
the liquid state is still the preferential one, but the time 
needed for a created fluctuation to vanish, increases dra- 
matically. At some temperature T(A), related to the 
factor A, we have tq(T) — > —00. We identify this tem- 
perature with the glass transition temperature T g , which 
depends on the cooling rate, as it should do. The con- 
figurational change r\ which may cause the relaxation for 
the low-temperature supercooled liquid has become in- 
finitely slow, thus the liquid behaves as a solid, and T g 
satisfies: 
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A(E C - E a ) 
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On the contrary, we can remark that starting from T = 
(at equilibrium) r (T) is close to r equ and the p(Co) = 
solution is now repulsive (fig. 2, solid line with circles). 
r (T) rapidly falls to 0, so the configurational change is 
possible and crystalline fluctuations can grow easily, in 
agreement with current observation stating that recrys- 
tallisation occurs when a glass is heated up. 
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Construction of the time-temperature-transformation 
(T-T-T) curves. As we have constructed this model by 
using a simple crystallization process, it is interesting to 
relate it to results concerning the kinetic treatment of 
glass formation, which describe the crystallization pro- 
cess by considering both nucleation and crystal growth 
and estimate the cooling rates required to form glasses 
P]. Thus, kinetic treatments of glass formation are based 
on identfying a certain value of the volume fraction crys- 
tallized v c /v, as borderline between an amorphous and 
a crystalline solid jll| (v c /v is generally of the order of 
the experimentally just detectable degree of crystallinity, 
i.e. around 10 -6 ). The results of such investigations are 
plotted on a t — T plane for different values of v c /v and 
represent the time-temperature-transformation (T-T-T) 
curves M- 




FIG. 3. The T-T-T curves for different s(t) behavior a) 
s(t) = 1 — e - *. Solid line (with points) To = 1.4, solid line 
To — 1.7, both with a fraction of x c = 1CF 6 . Dashed line: 
To = 1.7 and x c — 0.2. b) s(t) = t. Dots: system with 
To = 1.7, x c = 1(T 6 . As before, E c - E a = -1 for all situa- 
tions. 



In our mean-field description, the extensive variables 
p c and v c are related by: v c — p c v, where v is the total 
volume of the system. Therefore, if v c = p c v is main- 
tained constant, the right hand side of the dynamical 
equation will be equal to for any temperature. We can 
then compute the T-dependence of the factor A and ob- 
tain after integration a T — t relationship. If x c is the 
fraction of crystalline atoms (with probability p c nt ), then 



e 2a + e a 



A must satisfy: 
x c (e a - l) 2 
and integration yields: 



4<* 



a(x c - 1) 



a (x c - 1) - x c 



ote a dT 
~T~ s(t)dt 



= s(t) 



where T is the initial temperature of the melt and a the 
corresponding factor a. We have plotted different situa- 
tions in fig. 3 for different t-dependences of the function 
s(t). One natural (and the simplest) dependence for s(t) 
is the exponential one with s(0) = and s(oo) = 1, in 
order to agree with the construction of p c (t). Figure 3 
shows the dependence of the temperature with time if 
one keeps x c constant during nucleation. The curves dis- 
played show the same behavior as the constructed T-T-T 
curves by Uhlmann [|||l0| . i) The T-t curve with a lower 
degree of crystallinity (x c = 1CP 6 ) envelopes the curve 
with a higher degree [x c — 0.2), i.e. the temperature 
must decrease more rapidly, as seen on the figure, ii) 
The effect of the initial temperature To can also be ob- 
served. 

In conclusion, we should stress that the dynamical 
analysis presented here puts forward the general accepted 
picture that a glass is obtained by cooling a melt enough 
in order to avoid nucleation, although this picture dis- 
plays no transition in the thermodynamical sense. How- 
ever, different physical parameters are involved in this 
very complex transition such as free energy, cooling rate, 
nucleation frequency, initial temperature or structure. 
The work presented here was an attempt for the descrip- 
tion of the glass formation using a nucleation process. 
Forthcoming work will include structural factors pT[ . 



[1] W. Gotze, Science 267 (1995) 403 

[2] W. Bengtzelius, Phys. Rev A 34 (1986) 5059 

[3] Th.M. Nieuwenhuizen, Phys. Rev. Lett. 78 (1997) 3491, 

Phys. Rev. Lett. 79 (1997) 1317 
[4] F. Ritort, Phys. Rev. Lett. 75 (1995) 1190 
[5] J. Jackie, Phil. Mag. B 44 (1981) 533 
[6] D.R. Uhlmann, H. Yinnon, The formation of glasses, 

Glass, Science and Technology, Vol.1, Academic Press 

1983 

[7] see for example, T.R. Kirkpatrick, D. Thirnmalai; Phys. 

Rev. B 36 (1987) 5388 
[8] R. Kerner, D.M. Dos Santos, Phys. Rev. B 37 (1988) 

3881 

[9] R. Kerner, Physica B 215 (1995) 267 
[10] D.R Uhlmann, J. Non-Cryst. Solids 7 (1972) 337 
[11] M. Micoulaut, in preparation 



4 



